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ABSTRACT 

A technique is described for measuring electrical currents in the solar corona. 
It uses radioastronomical polarization measurements of a spatially-extended radio 
source viewed through the corona. The observations yield the difference in the 
Faraday rotation measure between two closely-spaced lines of sight through the 
corona, a measurement referred to as differential Faraday rotation. It is shown 
that the expression for differential Faraday rotation is proportional to the path 
integral § uB ■ ds where n is the plasma density and B is the coronal magnetic 
field. The integral is around a closed loop (Amperian Loop) in the corona. If the 
plasma density is assumed roughly constant, the differential Faraday rotation is 
proportional to the current within the loop, via Ampere's Law. The validity of 
the constant density approximation is discussed, and two test cases are presented 
in which the associated error in the inferred current is small, of order tens of 
percent or less. The method is illustrated with observations of the radio source 
3C228 with the Very Large Array (VLA) in August, 2003. A measurement of a 
differential Faraday rotation "event" on August 16, 2003, yields an estimate of 
2.5 X 10^ Amperes in the Amperian Loop. A smaller event on August 18 yields 
an enclosed current of 2.3 x 10® Amperes. The implications of these currents for 
coronal heating are briefly discussed. 

Subject headings: Sun: corona — Sun:magnetic fields — plasmas 



1. Introduction 

Electrical currents certainly flow in the solar corona. The structure seen in eclipse 
photographs or coronagraph images shows that pressure gradients must be balanced by elec- 
trodynamic forces. Although information on the strength and form of the coronal magnetic 
field is limited, it is clear that the true field is deformed from the potential field generated 
by a magnetic scalar potential. This difference is due to the presence of electrical currents. 
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In addition, there is a class of theories for heating of the solar corona which invoke Joule 
heating from c oronal currents , probably contained in turbulent current sheets. This idea 
originated with iParker I (Il972[). and has been elaborated in many subsequent works (see, for 
example , iGudiksen and Nordlund I l2005l : iPeter et al I l2006l ) . 



Observational measurements of coronal currents appear to be nonexistent. In fact, 
measurements of the coronal magnetic field itself at heliocentric distances of a few solar radii 
are limited to results from Fara day rotation of a trans-coro nal radio source and analysis of 
radio emission from solar flares (IBird and Edenhofer lllQQOl ). There is a literature on results 
from Faraday rotation observations of the large scale structure of the coronal field, as well 
as magnetic field inhomogeneities on a w ide range of sca l es. E x amples of such papers , whic h 
give references to the wider literatur e are iHoUweg et al I (Il982[): iBird and Edenhofer I (jl990l ) ; 
Mancuso and Spangler I Jl999l . boooh : ISpangler I J2005h and llnglebv et al I J2007h . 



In this paper, I discuss how Faraday Rotation measurements can also provide an obser- 
vational estimate of electrical cu rrents in the cor ona. The technique requires measurements 
of differential Faraday Rotation (ISpangler II2005I ). which is the difference in the Faraday ro- 
tation measure between two closely-spaced lines of sight through the corona. As discussed in 
the previously cited papers, a Faraday Rotation measurement yields the rotation measure, 
i?M, given by iKraus I (119661 ) 



RM 



UgB ■ dz 



LOS 



The fundamental physical constants of e,me,c, and eo are, respectively, the fundamental 
charge, the mass of an electron, the speed of light, and the permittivity of free space. The 
electron density in the plasma is ng, and B is the vector magnetic field. The incremental 
vector dz is a spatial increment along the line of sight, which is the path on which the 
radio waves propagate. Positive dz is in the direction from the source to the observer. The 
subscript LOS on the integral indicates an integral along the line of sight. Equation (1) is in 
SI units, as opposed to cgs, which has been used in our previous papers. The SI system is 
used in this paper for convenience in discussing electrical currents. The units of the rotation 
measure are radians/m^. 

As will be discussed in Section 2, when the rotation measures on two (or more) closely- 
spaced lines of sight are compared, one has an estimate of the electrical current between the 
two lines of sight. Such multiple lines of sight are available when one images an extended 
radio source (such as a radio galaxy or quasar) which is occulted by the corona. As is 
discussed in some of the papers referenced above, such measurements are straightforward 



- 3 - 



with the Very Large Array radiotelescope of the National Radio Astronomy Observator}0. 



The present investigation was motivated by laboratory experiments using differential 
Faraday rotation measurements to determine internal currents in the Madis on Symmetric 
Torus (MST) Reversed Field Pin ch (RFP) at the University of Wisco nsin ( iPrager I Il999l ) 



by D.L. Brower and W.X. Ding (IBrower et al I l2002l : iDing et al 1 120031 ) . In the case of the 



MST device, Ding and Brower use an infrared laser as the source of polarized radiation for 
measurement of Faraday Rotation. Measurement of the rotation measure along spatially- 
separated paths is achieved by directing the beam along distinct chords through the machine. 
In spite of obvious technical differences, the physical situation is virtually identical to that of 
the present paper. The inferred current profiles in the MST, and their dependence on time, 
are in agreement with theoretical pr e dictio ns ar id the resul t s from other diagnostics of the 
plasma. The results of iBrower et al I (120021 ) and iDing et al I (120031 ) encourage application of 
this technique to the solar corona. 

The organization of this paper is as follows. Section 2 describes the basis of the tech- 
nique, i.e. how differential Faraday Rotation measurements can provide a measurement of 
electrical currents in the corona. Section 3 provides an observational implementation of this 
technique with V ery Large Array observations of the radio source 3C228 in August, 2003 
(ISpangler 1 12005| ). In that section, it is shown that the total current contained with the 
Amperian Loop formed by the lines of sight was as high as 2.5 GigaAmperes during one 
period on August 16, but less than 0.8 GigaAmperes during a 3 hour period of high quality 
data before this event. A smaller, marginal detection of differential Faraday Rotation on 
August 18 yielded a current of the order of 0.23 GigaAmperes. In Section 4 we discuss 
the validity of an approximation used in the derivation of our expression for the current 
contained between two lines of sight through the corona. This assumption (stated in eq. [3] 
below) is that the plasma density can be approximated as roughly constant in the region that 
makes the maximum contribution to the Faraday rotation, and taken outside the integral 
defining the differential rotation measure. Arguments based on a simple theoretical model 
for coronal current sheets as well as data from a laboratory fusion plasma indicate that the 
error introduced by this approximation is relatively small, of order tens of percent or less. 
Section 5 briefly discusses the implications of our observations and analysis for theories of 
coronal heating by Joule dissipation of the currents. Our estimate of Joule heating using the 
Spitzer resistivity is about 6 orders of magnitude less than the observational estimate for 
the volumetric heating rate in the relevant part of the corona. Finally, Section 6 summarizes 



^The Very Large Array is an instrument of the National Radio Astronomy Observatory. The NRAO 
is a facihty of the National Science Foundation, operated under cooperative agreement with Associated 
Universities, Inc. 
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and concludes. 



2. Physical Basis of the Technique 



The technique is based on simultaneous Faraday rotation measurements along two lines 
of sight through the corona, which are separated by a small angular distance on the sky, and a 
corresponding physical separation / in the solar corona. Such observat ions can and have been 



made with the Very Large Array (VLA) of radio gal axies and quasars (ISakurai and Spangler 



I994J : iMancuso and Spangler Ill999l : ISpangler II2005I ) . A speci fic illustration is our observation 



of the radio galaxy 3C228 on August 16 a nd August 18, 20 03 (ISpangler II2005I : ISpangler et al 



20071 ). As may be seen in Figure 2 of ISpangler I (120051 ) . 3C228 is a double radio source 



with bright, highly polarized hot spots separated by about 46 arcseconds on the sky. The 
corresponding physical separation between the lines of sight to the two hot spots in the 
corona is about 33,000 km. The observations of August 16 and August 18 were made when 
the lines of sight passed within 6.7Rq and 5.2i?0 of the center of the Sun, respectively. 

A Faraday Rotation measurement of an extended radio source with two components is 
illustrated in cartoon form in Figure 1(a). The radio telescope measures rotation measure 
values RMa and RMb on the two paths which have a transverse separation /. The shaded 
area is meant to represent the coronal plasma, with the gray scale conveying the strength and 
sign of the current density. Black regions indicate regions of large positive current density, 
white areas are regions of large negative current density, and gray areas are regions of zero 
current density. Thi s picture is take n from numerical studies of current sheet development 
in MHD turbulence (ISpangler Ill999l ). 



The distribution of currents shown in Figure 1 is also consistent with a picture in which 
the current- carrying entities are coronal loops. This case would be relevant to streamer 
plasmas and the closed-field part of the corona. In this case, currents fiow along coronal 
l oops, up o n one footpoint and dow n on the other, as discussed, for example, by IChen 
( 119891 ) and I Chen and S chuck I (120071 ). The lines of sight to a radio source in a differential 
Faraday rotation experiment could "bracket" one branch of a loop and miss the other, so 
the cartoon in Figure 1 reasonably represents this well. 

The differential Faraday Rotation is defined as ARM = RM^ — RMb- As is illustrated 
in Figure 1(b), ARM is equivalent to the sum RMi + RM^, where RMi is the same as RMa, 
and RM3 is the same as RMb, but with the direction in which the integration is taken being 
reversed, i.e. from the telescope to the source. The value of ARM differs insignificantly 
from the sum of four terms, RMi + RM2 + RM^ + RM4, as shown in Figure 1(c). The 
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reason for this is that the contributions of terms RM2 and RM4 are small compared to that 
of RMi + RM3. There are two reasons for this. First, for the circumstances of a real coronal 
observation, the length of the segments 2 and 4 is small compared to 1 and 3. Second, these 
segments may be considered to be located far from the point of closest approach to the Sun 
(taken to be the middle of the diagram in all three panels), where both the plasma density 
and magnetic field strength are much smaller than their values near the point of closest 
approach. 

The net result then, is that the differential Faraday Rotation is given by 

ARM = RMi + RM2 + RM2 + RM3 jueB-ds (2) 
where the last expression represents a path integral around the closed Amperian Loop defined 

3 

by the segments 1,2,3, and 4 in Figure 1. The constant C is defined as C = STr^Jeom^ ~ 
2.631 X 10-^3 in sj nnii-s 

Equation (2) strongly recalls Ampere's Law, but the obvious difference is that in the 
present case the integrand in (2) is not just the magnetic field, but the product of the 
plasma density and the magnetic field. If we make the assumption that the measured ARM 
is dominated by a region in which the plasma density is relatively uniform, and given by a 
value n, we then have 

ARM = C jrieB-dsc^Cnj B- ds (3) 

and Ampere's Law can be utilized. The remainder of this paper will assume the convenience 
of equation (3). Although we think that equation (3) is a reasonable first-order approximation 
to employ, it is legitimate to ask for an estimate of the associated error, at least in some 
plausible test cases. This matter is discussed in Section 4 below. 

Given the approximation of equation (3), we can now use Ampere's Law, 

B-ds^ liol (4) 



where ^0 is the permeability of free space and I is the current contained within the Amperian 
Loop. Use of equation (4) in equation (3) allows us to write an equation giving the coronal 
current in terms of the measured ARM, 



Use of equation (5) requires an estimate of n. Fortunately, there are a number of 
empirical expressions for the coronal plasma density as a function of heliocentric distance. 
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For the purposes of this paper, we utihze the following expression ( ISpangler II2005I ). which 
has been employed in previous analyses of Faraday Rotation observations, and is in good 
agreement with independent estimates. 

\ -2.36 
112/ \ 



n{r) = 1.83 x 10'^ J m"^ (6) 

For purposes of simplifying the subsequent formulas, the index 2.36 will be rounded off to 
2.5. 

The density given by equation (6) would be the maximum that would be measured along 
a line of sight with a dimensionless impact parameter Rq = ^-^^ where r^m is the smallest 
heliocentric distance along the line of sight. The density n in equation (5) represents an 
average over the Amperian Loop defined by the observations. Therefore, for a given impact 
parameter Rq, n should be less that n{RQ) by some factor which is defined by n = an(i?o)- 
With these parameterizations, substitution of (6) into (5) yields the formula that will be 
used in Section 3, 

;(Ai?M)i?2-5 



'■obs 



1.65 X 10^^ '-^ Amperes (7) 



a 



The subscript "obs" on the current indicates that it is an estimate computed from observable 
quantities. Before leaving this section, it is worthwhile to emphasize the fundamental fact 
that the current given by equation (7) is the sum of all currents, positive and negative, within 
the Amperian Loop. It is thus obviously possible for strong currents to be flowing, but the 
total current given by equation (5) or (7) to be zero. 

The approach described above requires that the region which makes the dominant con- 
tribution to the Faraday rotation measure be sufficiently uniform to allow definition of a 
meaningful mean density n. As will be seen in Section 4.2 below, modest variations in 
plasma density do not result in a significant quantitative error when equations (5) or (7) are 
applied. However, larger errors could occur if one line of sight passes through a physically 
distinct region relative to the other. Whether this situation is pertinent will be determined 
by future observations which apply the technique described in this paper, and benefit from 
independent diagnostics of the plasma along the line of sight. 



3. Observational Implementation and Estimates of Coronal Currents 



We use data from the VLA observations of 3C228 in Augu st, 2 003. Further details o f 
the observations and data analysis are given in ISpangler I (120051 ) and ISpangler et al I (120071 ). 
The characteristics of the observations relevant to the present discussion are given at the 
beginning of Section 2 of the present paper. 
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Figure 3 of Spangler ( 2005 ) shows the Faraday Rotation % time series for 3C228 on 
August 16, 2003. The main feature of those observations was a large change in rotation 
measure over the 8 hour duration of the observing session. The data shown there distin- 



and a portion of the radio jet which hes between the hot spots (see Figure 2 of ISpangler 



2 of 


Spaneler 


(Spaneler et al 



20071 



These data have been used to calculate time series of ARM for both observing sessions. 
These data are shown in Figure 2, and are the differences in the rotation measures to the 
two hot spots of 3C228. The values of ARM plotted have used the data at both frequencies 
of observation, 1465 and 1665 MHz. In the case of the observations of August 18, 2003, only 
data from the last 4 hours of the observing session are used. Prior to this time, elevated 
system temperatures due to the proximity of the Sun caused a substantial loss in data 
quality. The most prominent sigr ial in the data fo r either day is the "event" at 19h UT on 
August 16. Th i s was discussed in ISpangler I (120051 ) and evidence for its credibility is given in 
Spangler et al I (120071 ). The differential Faraday Rotation had a maximum absolute value of 
6.52 ± 0.77 rad/m^ in the scan at 19h UT, and declined over the next two scans to a value 
around zero. 

The observations for August 18, when the line of sight passed closer to the Sun, show 
no such prominent event. The final scan shows a marginally significant detection of ARM = 
1.15 ± 0.42 rad/m^. Other measurements on this day are consistent with zero. 

Returning to the data for August 16, the three hours prior to the 19h event show 
measurements consistent with ARM = ± 2 radian/m^. These data can be used to set a 
significant upper limit to the current within the lines of sight. Data after 20h UT on August 
16 appear to show significant ARM values, although with no indication of a systematic 
mean value. We do not believe these data show good evidence for differential Faraday 
Rotation during this interval, because the errors are calculated from a propagation of known 
radiometer noise errors, and may (and probably do) underestimate the true errors due to 
solar interference. 

The final observational point to be made is that the ARM event at 19h on August 16 
occurred at a ti me when the rate of change of the rotation measure increased (see Figure 3 of 
Spangler II2005I ). As mentioned in that paper, the line of sight during this time was moving 
deeper into a coronal streamer, so this feature may not be typical of coronal properties. 



^The rotation in the polarization position angle, which is the quantity measured. It is given by the 
product of the rotation measure and the square of the observing wavelength. 



-8- 



In particular, it may reveal the system of currents associated with a streamer rather than 
a network of turbulent current sheets as illustrated in Figure 1, and which constitutes the 
model in the calculations to be briefly described in Section 5. 

The above data can be used in equation (7) to obtain estimates for the electrical current. 
The results are given in Table 1, in which results are given for 3 periods during August 16 
and 18, 2003. In the calculations used here, we have assumed a value of a (defined above, 
just prior to equation (7)) of i. In Table 2 we use the absolute value of ARM, since its sign, 
and the corresponding sign of /, are not considered in this paper. 
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Fig. 1. — Illustration of the measurement of coronal currents. Lines of sight through the 
corona are represented by dotted lines, with the arrow indicating the direction in which 
the integral is taken. The square shaded area represents the coronal plasma, with the 
gray scale indicating the coronal current density. Positive current density is black, and 
negative is white, (a) Illustrates the measurements of the Faraday rotation measure along 
two closely-spaced hues of sight, yielding RMa and RMb- (b) Illustrates the differential 
rotation measure Ai?M = RMa — RMb- It is equivalent to the sum of RMi and RM^. (c) 
The differential Faraday rotation measurement is very nearly the same as the sum of the 
four parts RMi + RM2 + RM^ + RM^, since the "end pieces" RM2 and RM^ contribute 
negligibly to the sum. The near-equality of the sum of the four segments and the differential 
Faraday Rotation measurement means that ARM is essentially equal to the path integral 
of nB around the Amperian Loop shown in panel ( c). Th e gray scale representation of the 
coronal plasma is taken from Figure 6 of ISpangler I (119991 ). 



20 

UT Time (hours) 



21.5 
UT Time (hours) 



Fig. 2. — Differential Faraday Rotation data for August 16, 2003 (left panel) and August 
18, 2003 (right panel). The measurements show the difference in the rotation measures to 
the hot spots of the radio source 3C228, which are separated by 46 arcseconds (33,000 km in 
the corona). The smaller number of measurements on August 18 results from reduced data 
quality prior to 20h UT, due to solar interference. 
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Prom Table 1 it can be seen that detectable differential Faraday rotation requires cur- 
rents of 10* — 10^ Amperes. Upper limits deduced from present radioastronomical observa- 
tions will then be in this range as well. 

Before leaving this section, we again emphasize that differential rotation measures ARM 
of even a few radians/m^ are rare in our two August, 2003 observing sessions spanning a 
total observing time of 16 hours. The apparent detection of coronal currents on August 
16 was an isolated event, and may be associated with a coronal streamer. The event on 
August 18 is of marginal statistical significance. Consequently, the typical magnitude of 
differential Faraday rotation through the corona at the heliocentric distances characteristic 
of these observations is probably less than that of the event on August 16, with net currents 
which are also correspondingly smaller. 

4. Validity of the Constant Density Approximation 

In this section we discuss the validity of the approximation employed in equation (3). 
The motiviation in adopting it is that it is a reasonable, first order approximation which 
permits application of Faraday's Law. The discussion in this section is organized into four 
topics. 

4.1. The Nature of Coronal Current Sheets 

It is obvious that until we have a better understanding of the physical nature of coronal 
current sheets, including their compressibility and structure, we will be unable to determine 
how close (3) is to an equality. The most obvious counterexample to (3), that in which the 
plasma density varies in a constant magnetic field, is not a satisfactory model for coronal 
plasma structures. In general, a spatial gradient in the plasma density results in a similarly 
directed gradient in the pressure. If the plasma structures are in approximate mechanical 



Table 1. Results on Coronal Currents 



Date 


UT Time (hours) 


ARM (rad/m^) 


/ (Amperes) 


Aug. 16 


15.5 - 18.5 


< 2.0 


< 7.7 X 10* 


Aug. 16 


19.0 


6.52 ±0.77 


(2.50 ± 0.30) X 10*^ 


Aug. 18 


22.7 


1.15 ±0.42 


(2.34 ± 0.86) X 10* 
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equilibrium, this pressure gradient must be balanced, and the most obvious possibility is a 
J X B force, where J is the current density. Thus currents will be necessary where there are 
plasma density variations, and while the two sides of equation (3) might not be equal, they 
should be of comparable magnitude. 

A case in which the above considerations are not valid would be if differential Faraday 
rotation signals are due to parallel-propagating slow magnetosonic waves, for which there are 
density variations but no field variations. There is also the possibility that density variations 
are associated with anticorrelated temperature variations, so that pressure variations are 
absent, and the need for currents to achieve mechanical equilibrium is removed. 



4.2. The Z-Pinch as a Model for Coronal Current Filaments 



In this subsection we adopt a simplified model for coronal current filaments to explicitly 
calcula te the error associated with the approximation (3). The mod el we adopt is the Z- 
Pinch (IKrall and Trivelpiece Ill973l : iGurnett and Bhattacharjee II2005I ) which is a solution to 
the equation for mechanical equilibrium 

JxB-Vp = (8) 

in cylindrical symmetry, where p is the plasma pressure. In the Z-Pinch, current can flow 
in the z direction (coordinate perpendicular to the plane of the page in Figure 1, not the 
coordinate of the observational line of sight defined in eq.[l]) and 6 direction (azimuthal 
direction in the plane of Figure 1), and thus balance pressure gradients in the radial direction 
in this plane. The Z-pinch thus resembles the current sheets shown in Figure 1, with the 
exception that we now approximate these current filaments as having circular symmetry. 



Although the Z-Pinch is an unstable equilibrium (IKrall and Trivelpiece Ill973l : iGurnett and Bhattachar 



20051 ). it is a well-studied model for current-carrying plasma structures, and one can imagine 
compressible astrophysical turbulence as consisting of such structures forming, persisting for 
a while, then dispersing with subsequent reformation. 

We adopt the common simplification of ignoring the 9 component of the current, which 
is equivalent to assuming no radial gradients in the z component of the magnetic field. 
Equation (8) then yields the following differential equation for the relation between the 
plasma pressure and the 9 component of the magnetic field 



Be djrBe) 
r dr 



dp 
dr 



(9) 



where r in this section is the radial coordinate in a cylindrical coordinate system, not 
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the heliocentric distance as used in Section 2. It is Bg that causes the differential Fara- 



day rotation in this model . In discussions of the Z-Pinch (IKrall and Trivelpiece I Il973 



Gurnett and Bhattacharjee II2005I ). one prescribes one of the fields p or Bg, then uses equa- 
tion (9) to determine the other. This procedure is well-suited to our present task, which 
is to see how spatial gradients in the plasma density of a plausible plasma structure affect 
differential Faraday rotation. 

We assume an isothermal approximation for the plasma equation of state, which will 
give a larger density change for a given pressure change than adiabatic equations of state. 

p{r) = c^p = clmnn^r) (10) 

where Cs is the sound speed, p is the mass density, mn is the mass of a hydrogen atom, 
and n(r) is the number density of atoms. We assume a pure hydrogen plasma. Substituting 
equation (10) into (9), and making the change of variable Bg —> rBg = /, equation (9) 
becomes 

f!W ^ _27r^^ (11) 
dr dr 

with 7 = fiQC^niH, which is directly integrable. The solution is 

r dn 

fir) = {rBe{r)f = -2^ j r^—dr (12) 

We adopt a model for the density profile which is enhanced in the current channel, 

n(r) = n + me-^"/") (13) 

with a being the specified radial extent of the current filament. We substitute equation (13) 
into (12), and change variable of integration from r ^ ^ = ^. The solution then becomes 

f\r) = {rBg{r)f = 2^nia^F{0 (14) 

F(0 = / ee~^d^ = 2-{e + 2^ + 2)e^« (15) 

Be{r) = ^^a^/F(0 (16) 
r 

where Ap = c^mHUi. Equation (16) also yields the expression for the current density 

Mr) = (^) ^/^G(0 (17) 
G{0 - , 
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We now express the different fields in tfie current filament in terms of common param- 
eters which characterize the host plasma and the current channel. The plasma pressure is, 
from above, 

p — clniHin + Tiie~^) (19) 
We define the "modulation index" e = ^ = so that 

n^n{l + ee-^) (20) 
p = po{l + ee-«) (21) 

With Po = clmnn. We now relate the plasma pressure to the magnetic pressure via 

PO = P{^) (22) 

Ap = el3A (23) 

Where (3 is the conventional plasma (3, here defined as the ratio of gas pressure to magnetic 
pressure, and Bq is the mean magnetic field. Using equation (23) in the equation for Bg 
gives 



Be = BoVePF2{0 (24) 
F,{0 = (25) 

With these expressions, we can now calculate both sides of the near-equality (3). The 
integral is taken in a clockwise sense, as indicated in Figure 1(b), on the Amperian Loop 
shown in Figure 1, with the convention that the origin of the polar coordinate system is 
midway between the integration paths A and B. The results of the calculation below are not 
qualitatively changed if the center of the current filament (r = 0) is not midway between the 
two linear paths of integration A and B. The segments of the Amperian Loop connecting 
paths A and B (segments 2 and 4 in Figure 1(c)) are assumed to lie at infinity and are 
neglected in the calculation below. 

The integral around the Amperian Loop is then given by symmetry considerations as 

nB-ds = 4 I nB- ds (26) 



6 



where the path segment b corresponds to | > ^ > 0. In what follows we evaluate the integral 
over path b. Given the expression for Bq, etc developed above, the integrand of equation (3) 
is 

nB-ds^ n{$)B0{^)e0 ■ e^dx (27) 
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where is a unit vector in the the x direction (horizontal axis in Figure 1), and eg is the unit 
vector in the direction of increasing d. The dot product of the unit vectors S&e^-ex — — cos ^, 
where Q is the complement of Q. 

An important variable is the perpendicular distance it! of the integration path from the 
center of the current channel. In terms of our observations, 2R would be the separation / 

between the hnes of sight through the corona. The incremental path element dx = R sec^ 6d6, 
and the variable ^ ~ ^ sec 9. Substituting all of these relations into (27) and forming 
the integral, we have 

^nB-ds^Aj^ n{C}Be{i)Rsecede (28) 
Substituting equations (20) and (24) into (28), we have 

(29) 



nB -ds^AnoBo^R / Fa (0 sec + e / e-^F2{C) secM 

Jo Jo 



^AnoBoVePR[Sa + eSh\ (30) 



= AnoBo^/€(3RSa 



(31) 



Equations (29) and (30) define the integrals Sa (f ) and Sa (f ). An appealing feature of 
equation (31) is that the leading multiplicative term on the right hand side is the value of 
f nB • ds for the case of uniform density n. The second term in brackets in (31) is therefore 
the correction due to the change in density in the Z-Pinch. It is directly proportional to the 
density modulation index e. 

The parameter F = ^ gives the fractional error in approximating f nB-ds hy n f B-ds. 
Obviously, this is equal to the fractional error in the current inferred from the radioastro- 
nomical technique which is the subject of this paper. Given the definition of the density 
modulation index e, one would expect it to be less than unity, say e < 0.50, although this 
is not strictly required in the definition of equation (20). In any case, it is a parameter 
which might be observationally accessible. Given an adopted value of e, the correction factor 
which relates the true current in a plasma to that retrieved by the technique of Section 2 is 
then dependent on the ratio The integrals Sa and Sj, were evaluated numerically with a 
Mathematica notebook. The results for a range of values of ^ are given in Table 2. 

The conclusions to be drawn from the calculations summarized in Table 2 are clear. 
Even for hnes of sight which pass deep in the interior of the current channel, the ratio 
1^ < 1 by a substantial factor. The fractional error F = ^ will be smaller still, even for 
relatively large compression factors of e ~ 0.5. For example, if the line of sight passes within 
i? = 0.1 — 0.2a of the center of the filament, and the compression factor e = 0.50, the error 
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in the current retrieved by the method of Sections 2 and 3 would be only 17 - 19 %. This is 
completely negligible in the context of this paper, and indicates that the method described 
here would be quite accurate. 

For lines of sight that pass further from the center of the filament, the fractional error 
becomes smaller still, approaching zero as — ^ 1, as indeed it must. 



4.3. Analysis of Measurements from the MST 



In this subsection, we consider a test of our method with data from t he Faraday rotation 



experiment on th e Madison Symmetric Torus (MST) plasma machine ( iBrower et al 1 12002 



experiment 


on tl 


Dine et al 


2003) 



contrast to the mathematical ideal of Section 4.2. In addition, as a laboratory experiment, 
it is diagnosed in a way that the solar corona cannot be. 

In some sense, use of the MST data assures corroboration of the technique. The goal of 
the experiments was to measure currents flowing in the MST, and those ex peri ments moti- 



vated the current investigation. However, the analysis of iBrower et al I (|2002| ) and lDing et al 



( 120031 ) utilized more experimental information than is available to us, and (at least in part) 
applied more sophisticated methods of processing the data, which benefitted from a-priori 
information on the structure of the toroidal plasma. 

In this subsection, we app ly the same formul as used in Sections 2 and 3 for 3C228 to 



published data from the MST ( Brower et al 20o3). and com pare the retrieved current with 



that known to flow in the MST device. IBrower et al I (120021 ) report that a current of 400 kA 
flows in the MST during their experiments, and that the mean density (which we adopt for 
n) is 10^^ m~^. The wavelength of the polarized laser beam used to probe the plasma was 
432 nm. 



We utilize the data shown in Figure 1 of IBrower et al I (|2002| ). which shows Faraday 
rotation of about 5° for the chord at an offset distance of -17 cm (corresponding to the 
parameter R defined in the previous subsection), and Faraday rotation of about —3°. 5 for 
the chord at i? = +21 cm. We thus have differential Faraday rotation of Ax ~ 8°. 5, 
corresponding to a differential rotation measure ARM ^ 8.0 x 10^ rad/m^ between two lines 
of sight separated by 0.38 meters. Substitution of this value for ARM into equation (5) 
yields a current of 241 kA. 



This value itself, resulting from the most straightforward application of the methods 
used in the study of the corona, is in reasonably good agreement with the known value. 
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It is within a factor of two of the known, measured current in the MST device. Similar 
precision in the measurement of coronal currents would be cause for delirium. However, we 
can improve agreement between remotely inferred currents in the MST, using equation (5) 
of this paper, and the known value for the current. 



The laser polarimeter system of iBrower et al I (120021 ) and iDing et al I (120031 ) utilized a 
number of chords with different values of R, so the current profile can be retrieved. In fact, 
the change in the f unctional form of the current density during a sawtooth crash on MST 



was a main result of Brower et al 



Brower et al 



(120021 ) . The inferred toroidal current density profile Jz(r) 



(I2OO2I ). and shows that the current channel has a 



is plotted in Figure 4 of 
greater radial extent than 0.2m. Accordingly, the estimate of I immediately above must 



be an underestimate. T o approximate 



y cal culate the current within the two outer chords 



displayed in Figure 1 of iBrower et al I (120021 ). we represent the current density profile Jz{r) 
by a Gaussian, 



JJr) 



J(0)e 



(32) 



We use plot 4(b) of IBrower et al I (120021 ) and the profile measured before a sawtooth crash, 
for which we estimate values of J(0) = 1.93 x 10^ A/m^, and a = 0.29 m. The total current 
within two lines of sight separated by a distance 2R with the center of the current filament 
midway between the lines of sight is 



I{R) = TTj{0)a^Erf{X) 



(33) 



where Erf{X) is the Error Function oi X = ^. Using the above values of J(0) and a, 
and with R = 0.19 m, equation (33) gives I{R = 0.19m) = 331kA. This value is, as 
expected, in better agreement with the value inferred from equation (5), and indicates a 37 
% error between the "remote sensing" and "true" value, a difference which includes possible 
systematic errors in the plasma physics experiment, which are obviously beyond the scope 
of this paper. In any case, the degree of agreement supports the credibility of the technique 
discussed in this paper. 

The work of Brower and Ding demonstrates that future coronal differential Faraday 
rotation observations, with measurable ARM on several lines of sight to a background radio 
source, could measure not only the total current but the profile of coronal current density, 
as is done in the MST. To summarize the results of this and the preceding subsection, 
investigation of a simplified analytical model of coronal current filaments and experimental 
results from a laboratory plasma have suggested that the approximation of equation (3) 
results in errors in the retrieved current of, at most, tens of percent. Such an error does not 
detract from the significance of our results for coronal physics. 
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4.4. Prospects for Improved Modeling of Coronal Currents 

The analysis of Section 4.2 would be improved through the use of realistic models for 
the compressive current systems which exist in the solar corona. As in the case of Sections 
4.2 and 4.3, the goal will be to determine the error associated with assuming a constant 
density in the plasma, so that the associated error in the retrieved current can be estimated. 

A promising set of exis t ing ca lculat ions is comprised of the recent coronal simulations 



of iGudiksen and Nordlund I (120051 ) and iPeter et al I (120061 ) , which are specifically directed 



to studies of Joule heating of the corona. The calculations presented in those papers are 
three dimensional and give the vector magnetic field, plasma density, and current density at 
all locations in the simulation. The Faraday rotation measure along different lines of sight 
through the simulation could be calculated and the relationship between differential Faraday 
rotation and current properties investigated. 

The investigations to date by Gudiksen, Peter, and collaborators have concentrated on 
the corona at much smaller heliocentric distances than probed by our VLA observations. 
However, future calculations could presumably be carried out to simulate conditions at 
heliocentric distances of 5 — IORq. Furthermore, even the present calculations which have 
been described in the literature could be used to examine the general relation between 
differential Faraday rotation and current filaments in a simulated environment patterned 
after the corona. 



5. Implications for Coronal Heating 

Having provided observational estimates of coronal currents in Section 3 of this paper, 
there arises the obvious question as to whether these currents can provide the requisite 
heating of the corona via Joule heating. Unfortunately, it is not straightforward to answer 
this question. The volumetric heating rate, which must be the basis of such a discussion, 
is r/J^ where r] is the electrical resistivity of the plasma, and J is the current density. To 
obtain the average coronal heating rate in the region probed by our observations, one must 
estimate the geometry, extent, and number of the current sheets in a given volume of the 
corona, in order to estimate both J from a measurement of J, and also to calculate the filling 
factor of the volume in which heating is going on. This obviously requires a model for the 
current sheets. Furthermore, and most importantly, an expression for the resistivity r] must 
be chosen. 

We carried out a calculation for the average volumetric heating rate, using two models 
for the nature of the current sheets responsible for the differential Faraday rotation. These 
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calculations envisioned the corona as illustrated in Figure 1, with thin, long (along the large- 
scale coronal field) sheets, in which the sign of the current density would vary from one sheet 
to the next. One of the two models assumed that there was equal probability of positive and 
negative current density, as would be expected in a true turbulent situation, and in which 
the measured differential Faraday rotation is due to a statistical fluctuation in the number 
of positive and negative sheets within the Amperian Loop. The second model assumed a 
statistical preference for one sign of the current density. 

Formulas were derived for effective volumetric heating rates in both cases. These heating 
rates depend on assumed properties of the current sheets, such as thickness, lateral extent, 
extension along the large scale fleld, etc. The heating rates also depend crucially on the resis- 
tivity. Plausible estimates could be obtained for all the g eometric properties of these current 



sheets. For the resistivity, we used the Spitzer resistivity (IGurnett and Bhattacharjee II2005I ) 
which is based on fully-understood principles and may be assumed to be a lower limit to the 
true coronal resistivity. The volumetric heating rate for either mode l was r oughly 6 orders 



of magnitude less than that estimated by lCranmer and Ballegooijen I (120051 ) as necessary to 



explain heating in the relevant part of the corona. This calculation indicates that either the 
currents which are reported here are irrelevant for coronal heating, or that they do play a 
role in coronal heating, but that the effective resistivity of the coronal plasma at heliocentric 
distances of 5 — IORq is approximately 6 orders of magnitude greater than the Spitzer value. 
A document with a detailed description of these calculations is available from the author. 



6. Summary and Conclusions 



We have pointed out that a type of polarization measurement (differential Faraday 
Rotation measurements on an extended radio source) which can be done with the 
Very Large Array can yield estimates of electrical currents in the solar corona. The 
technique could, in principle, be done with some other radio telescopes as well. This 
technique is an astronomical ada ptati on of a dia g nostic used in fusion plasmas, and 
described by iBrower et al I (120021 ) and iDing et al I (120031 ) . 



2. The technique has been applied to observations of the radio source 3C228 which were 
made on two days in August 2003, when the radio source was viewed through the corona 
at "impact parameters" of 6.7 and 5.2Rq. Detectable differential Faraday rotation 
was detected in relatively brief "events" on both days, yielding estimated currents of 
2.5 X 10^ and 2.3 x 10^ Amperes, respectively. Another interval of high quality data 
on one of the days yielded an upper limit to the differential Faraday rotation, and a 
corresponding upper limit to the current of 8 x 10*^ Amperes. 
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The technique as defined in this paper is dependent on an approximation in which the 
plasma density is assumed constant in the region within the Amperian Loop which 
makes the dominant contribution to the rotation measure. The vahdity of this approx- 
imation was investigated in two ways in Section 4. First, the value of f nB ■ ds was 
compared with n f B-ds for a Z-Pinch model in which the plasma density is modified in 
the current filament. Second, the formula used in our radioastronomical observations 
(eq.[5]) was applied to measurements made with an infrare d laser polarimeter on the 



MST plasma machine (IBrower et al II2002I : iDing et al 1120031 ) . and used to estimate the 
current, which could also be measured independently. The values agreed to 37 %. 

4. We explored the implications of the measured currents for Joule heating of the corona. 
The magnitude of the estimated, average volumetric heating rate is rendered uncertain 
by poorly-known properties of the current sheets or filaments, and by the value of the 
coronal resistivity. Use of the Spitzer resistivity leads to a value for the volumetric 
heating rate approximately 6 orders of magnitude less than the value indicated by 
observations. We conclude that either these currents are irrelevant for coronal heating, 
or that the true resistivity in the corona exceeds the Spitzer value by several orders of 
magnitude. Resolution of this matter obviously lies in a better understanding of the 
resistivity in a collisionless plasma. 



This work was supported at the University of Iowa by grants ATM03-11825, and ATM03- 
54782 from the National Science Foundation. I thank W.X. Ding and D.L. Brower of the 
University of California, Los Angeles, for interesting and helpful discussion about Faraday 
rotation diagnostics of the current profile on the MST Reversed Field Pinch device of the 
University of Wisconsin. 
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Table 2. Validity of Constant Density Approximation 
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Sa 




Sb 

Sa 
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0.645 
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1.585 


0.542 


0.342 


0.4 


1.443 


0.397 


0.275 


0.6 


1.331 


0.297 


0.223 


0.8 


1.237 


0.224 


0.181 


1.0 


1.154 
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1.015 
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0.902 


0.057 
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520 Edgemont Road 



C/^ ; Charlottesville, Virgina 22903 



Dear Tim 



I have just submitted to the Astrophysical Journal server the revised version of my paper, "A 



. Technique for Measuring Electrical Currents in the Solar Corona", manuscript ^ 71459. To 



review the history of this paper, it was submitted to the Astrophysical Journal on February 



O ! 16 of this year. I received the report of the referee (which included comments both by the 



referee and an observational colleague) on February 26. Your assessment of the referee's 
report was that the paper was "unacceptable in its current form". 



I was encouraged to revise and resubmit this paper for two reasons. First, the referee was 
2 ! favorably inclined to the basic idea of the paper, saying that "the idea is explained clearly, 
^ I and the paper should be published in some form" , and "the central idea of this paper should 
be published". Second, it seemed clear to me that the referees' objections were mainly in 
two areas which could be addressed. It is my hope that I have done so. 



The revised paper is a major revision, and involves a substantial amount of new work, under- 
taken specifically to address the criticisms of the referees. To guide the referees in assessing 
the revised version, and as an overview to what I have done, the first three sections (up 
to and including the calculation of currents from the observations) are only slightly revised 
from the original, in ways which incorporate suggestions from the referees. A completely new 
Section 4 has been inserted to address the validity of the constant density approximation. 
Finally, the revised Section 5 is, as recommended by both referees, a "severely abbreviated" 
version of the calculation of the heating rate by Joule heating of these currents. 

The criticisms of the primary referee are contained in the "Three Disappointments" , a term 
that suggests a dark variation on an entree in a Chinese restaurant. I hope my responses 
to the "Three Disappointments" fully and positively respond to these suggestions. There 
are also the criticisms and suggestions by the observational colleague. The observational 
colleague's remarks consist of a preamble, three enumerated comments on the method, and 
two comments on the interpretation. 



In what follows, I give my response to all of these points. 



The "Three Disappointments" 



(1) The referee wishes to have a discussion of the approximation/assumption of pulhng 
the density outside the path integral in equation (3) (both versions). Although I thought 
this was a reasonable first order approximation, the referee obviously has a legitimate issue 
here. In response to this point, I have included an entire new section, the present Section 
4 of the paper, entitled "Validity of the Constant Density Approximation" . There are four 
points brought up in this section. (A) Until wc have a better idea of the structure and 
thermodynamics of coronal irregularities, we cannot really answer this question. However, I 
do make the point that if there are density fluctuations, there will generally (but not in all 
cases) be corresponding electrical currents. (B) The main part of this section. Section 4.2, 
is a quantitative discussion of this approximation for a model of coronal current fllaments 
based on the Z-Pinch. The Z-Pinch is a current- carrying filament which confines an enhanced 
density hy J x B forces. This model is analysed quantitatively and a parameter is obtained 
(the parameter F introduced in the last paragraph on pl5) which gives the error due to the 
approximation involved in equation (3) of this paper. The quantitative results, given in Table 
2 of the revised version, show that the error should be of order 20 % in a fairly unfavorable 
case, and much less otherwise. Obviously, we have no assurance that the coronal current 
sheets are structured like Z-Pinches, but the point here is that in the case of a plausible 
facsimile for these loops, the error in the constant density approximation is very small. (C) 
I took data from the figures in the MST Reversed Field Pinch experiments (the experiments 
which inspired the development of this technique) and applied my formulas which were used 
for the observations of radio galaxies through the corona. This gave a "remote sensing" 
estimate of the current in the machine, which is independently known. In Section 4.3 of 
the revised paper, I describe these calculations, and show that the "retrieved" current is 
within 37 % of the known current in the plasma device. So this technique is shown to work 
(at the level of a few tens of percent) in the case of a "real" plasma. (D) I point out that 
this issue could be settled definitively by calculating the differential Faraday rotation along 
paths through numerical simulations of the corona in which Joule heating is playing a major 
role, and I reference the papers of Gudiksen, Nordlund, and Peter. This would be a very 
straightforward thing to do, and I hope my paper (if it is published) could encourage such a 
project. 

(2) The referee asked that the possibihty be mentioned that the currents could be those 
directed along coronal loops rather than the set of turbulent current sheets that I had 

envisioned. In the revised paper. I have included a paragraph of new text, third complete 
paragraph, p4, in which I discuss this possibility and reference two papers of Jim Chen's 
which specifically discuss (and illustrate in the figures) electrical currents fiowing on coronal 
loops. 

(3) The referee obviously did not like the section on estimates of Joule heating at a rate 
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determined by the Spitzer resistivity, and stated in his or her recommendation # 4 that 
"the discussion of coronal heating by the putative current should be severely abbreviated 
or deleted entirely". 1 was very reluctant to do this. 1 anticipated at the outset that Joule 
heating at the Spitzer rate would be far smaller than real rates, but 1 felt it was important to 
actually demonstrate that this is the case. Furthermore, although the final expression for the 
heating rate involved many "imponderables" , I think it important to go through the exercise 
and identify what parameters are most important in determining the true heating rate. The 
referee may be interested in knowing that I undertook this calculation at the encouragement 
of Peter Gary. 

Nonetheless, Peter Gary was not the referee of this paper, so I felt I had to modify my paper 
to satisfy the current referee. 1 have therefore chosen the "severely abbreviated" route, and 
the new Section 5 consists of only three paragraphs of text and no equations, and summarizes 
the results of the old Section 4 of the original paper. I also mention that the details of this 
calculation are available from the author. I really think I have to do this. In view of the 
prominence of Joule heating theories for coronal heating, a report which describes how some 
observational handle can be had on currents has to say something about how these numbers 
relate to Joule heating calculations. 

Comments of the Observational Colleague 

The preamble to this person's remarks makes the entirely valid point that the main "event" 
on August 16, 2003 occurred at a time when the line of sight to the radio source was passing 
into a streamer, so the currents which are being detected might be rare and associated with 
a single structure, rather than due to turbulent current sheets of the sort I have envisioned. 

I agree with everything the referee says here, and believe I clearly described the situation 
in the original version of the paper. However, to emphasize this point, I have inserted a 
paragraph of new text on pll, second complete paragraph, last paragraph of Section 3, 
which clearly indicates that most of the time, the differential rotation measure was less than 
that observed at 19h on August 16. In this I indicate, the sort of analysis here could 

be used to set upper hmits to the current. With regard to the marginal nature of the ARM 
on August 18, I stated the error on the measurement in the original version as well as the 
current version, so I did "put my cards on the table" . Nonetheless, to make this completely 
clear 1 have appended the modifier "marginal" when referring to it in the revised version of 
the paper (see, for example, the description of this measurement in the last paragraph of 
Section 3). Once again, one may certainly use this measurement for an upper limit to the 
current. 

(4) The first comment in "issues with the method" notes that the differential Faraday 
rotation measurement is an integral over the path through the corona. I would respond that 
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this is completely obvious from the equations that describe the observations, beginning with 
equation (1). Furthermore, Ampere's Law actually requires that the measurement be related 
to an integral around the whole curve (the Amperian Loop). The point that one is averaging 
currents that are both small and large, and of both signs does need to be kept in mind, as 
is pointed out in the two sentences following equation (7), revised version. 

(5) The referee makes the very good point that if one line of sight goes through a region 
which is physically distinct from that perforated by the other LOS, a differential Faraday 
rotation signature would be observed which would be hard to interpret. I convey this point 
in a paragraph of new text inserted on p6 at the very end of Section 2. I also point out that 
future observations, using radio data such as mine and independent coronal diagnostics such 
as coronagraph observations, might indicate if this extreme inhomogeneity is affecting the 
observations. However, as 1 also point out in this new paragraph, the analysis of Section 4.2 
indicates that mild density inhomogeneity does not invalidate the estimates of the current. 

(6) I don't agree with the assertion that this method is reliant on assumptions about lack 
of time variability. Ampere's Law does not contain time derivatives (assuming that the 
displacement current is negligible), and the VLA makes observations along the two lines of 
sight A and B simultaneously, i.e. no time lag at all. We do not rely on Taylor's hypothesis 
of frozen plasma fluctuations as is sometimes done in analysis of spacecraft data. 

(7) Comments A and B, as well as their summary at the end of the report raise questions 
about the plausibility of parameters of the current sheets in my Joule heating calculations. 
The referee speculates whether the heating rate is based on a large number of contrived 
parameters. 1 can give a defense of the estimates 1 chose, as well as a frank "1 don't know" 
when it is warranted. However, since the section containing these calculations has been all 
but eliminated, on the basis of the recommendations of both referees, I think this point is 
now moot. 

The new, extensive discussion of the validity of the constant density approximation, as well 
as the de-emphasis of the coronal heating calculation, are reflected in the revised abstract 
and conclusions. 

This concludes my response to the referees. I hope that it is clear that I have taken their re- 
marks seriously, and substantially modified the paper in response to their recommendations. 
I have done this even though I would have preferred to retain the calculation on heating. I 
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hope that the paper is now acceptable for pubhcation in The Astrophysical Journal. 

Sincerely yours, 
Steve 



Steven R. Spangler 
Professor of Physics and 
Astronomy 
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